Effective potentials for atom-atom interaction at low temperatures 
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We discuss the concept and design of effective atom-atom potentials that accurately describe 
any physical processes involving only states around the threshold. The existence of such potentials 
gives hope to a quantitative, and systematic, understanding of quantum few-atom and quantum 
many-atom systems at relatively low temperatures. 

PACS numbers: 05.30.-d,34.10.+x,03.75.Fi,34.20.-b 



The concept of model potential has played an im- 
portant role in many branches of physics. The well 
known examples include the Morse potential jjj and the 
Lennard- Jones potential for molecular systems, and the 
hard-sphere potential and the delta-function pseudopo- 
tential widely used in many-body theories, including 
theories for Bose- Einstein condensates (BEC) ra. 

There are many reasons why one uses a model potential 
instead of the "real" potential. Here, we only empha- 
size that a problem can simply become unmanageable 
if the "real" potential is used. For quantum systems, 
this statement quickly becomes true for three or more 
atoms. This explains, in part, that despite its many lim- 
itations, it has proven difficult to go substantially beyond 
the Gross-Pitaveskii theory for BEC that is based on the 
delta-function pseudopotential || . 

Our goal here is to discuss the concept and design of 
model potentials that better reflect the reality of atom- 
atom interaction than either the hard-sphere potential or 
the delta- function pseudopotential, yet simple enough to 
allow for applications in quantum few-atom and quantum 
many-atom systems. Here, being "simple" , to a large 
extent, means being shallow, as it is the depth of a po- 
tential, which can be measured by the number of bound 
states it supports, that determines the complexity of the 
resulting quantum few-atom and quantum many-atom 
problems. 

One of our key conclusions is the following. For N- 
atom states around the N-atom threshold (such as the 
BEC state [||), or physical processes that involve only 
states around the threshold (such as the three-body re- 
combination process ||), the interaction potential be- 
tween a pair of atoms, no matter how deep it might be, 
can be replaced by an effective potential supporting only 
one or a few bound states. Furthermore, because dif- 
ferent partial waves are described by the same effective 
potential, it can be used in precisely the same manner 
as any "real" potential. And in doing so, one reduces 
the complexity of the resulting quantum few-atom and 
quantum many-atom problems to a level comparable to 
that for He, a level that we are quickly learning to deal 
with §,§. 



If one thinks of the delta-function pseudopotential 
Q as describing the atomic interaction at the longest 
length scale in the zero-energy limit, 2-7r/fc, the nat- 
ural next step is the description of atomic interac- 
tion at the next, shorter, length scale. This scale is 
j3 n = (2/iC n /Ti 2 ) 1 ^ n ~ 2 \ which characterizes the long- 
range atomic interaction of the form of —C n /r n (n > 2). 
The angular-momentum-insensitive quantum-defect the- 
ory (AQDT) [Q, D provides a systematic description of 
atom-atom interaction at this scale, and is the basis of 
our discussion. 

Reference fffl focused on two-atom systems with 
V(r) — > —Ce/r e , but the same concepts and formula- 
tion are readily generalized to any n > 2. In this general 
formulation, a two-atom system with an asymptotic po- 
tential of the form of —C n /r n (n > 2) is described by 
a dimcnsionlcss K matrix K c (e, I) and a set of universal 
functions that are determined from the solutions of 



dr 2 



Ki + i) , K 



+ 1 



Uel{r) = 0, 



(1) 



where e = 2fj,e/h 2 . Specifically, K c (e,l) is defined by 
writing the wave function at large distances as a linear 
superposition of a pair of reference solutions of Eq. (pi): 



u el (r)=A el [fc l (r)-K c (e ) l)g c d (r)} 



(2) 



where f c and g c are purposely chosen to have the behav- 
ior 

f d ^ (2A) 1 /2( r//3 „)«/4 cos(2/ _ 7r/4) ; (3) 

ff c r^ n _ (2/7r) l/2 (r//3n) n/4 sin {y _ ^ (4) 

for all energies @. Here y = [2/(n - 2)] (/?„/r) ( "~ 2)/2 - 

AQDT asserts that K c (e, I) is approximately a con- 
stant that is independent of both e and /, provided 
that p n is greater than other, energy- independent, length 
scales present in the system Q]. For our purposes here, 
the most important conclusion of AQDT is the follow- 
ing. To the extent of K c (e, I) being energy and angular- 
momentum independent, potentials with the same type of 
long-range behavior (namely the same n) and the same 
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K° have, on a scaled energy basis, the same bound spectra 
and scattering properties around the threshold [Q. Here 
the bound spectra and scattering properties include all 
angular momentum states for which K c remains approx- 
imately ^-independent. 

Ignoring scaling relations implied by this statement, 
which will discussed elsewhere ||, AQDT gives the fol- 
lowing simple prescription for designing an effective po- 
tential that has the same physical properties around the 
threshold as the system of interest. First, choose a model 
potential. The only restriction is that it should have 
the right asymptotic behavior. Second, adjust the short 
range parameters of the model potential so that 



K c mod (0, l) = K c (0,l), 



(5) 



for one particular I. These two conditions do not uniquely 
determine an effective potential. Another auxiliary con- 
dition, which gives a convenient characterization of the 
depth of a potential, is the number of bound levels sup- 
ported by a model potential for a particular I, JVj. For 
classes of model potentials discussed below, these condi- 
tions uniquely determine a model potential. 

We stress here two classes of model potentials for which 
K c (0, 1), and the number of bound levels for each I, Ni, 
can be found analytically. One class is of the type of a 
hard-sphere with an attractive tail: 



Vhst{t) 



oo , r < tq 

-C n /r n , r > r 



(0) 



which will be denoted by HST. The other class is of the 
type of Lennard- Jones LJ(n, 2n — 2): 



V LJn (r) = -C n /r n + C 2n - 2 /r 



2n-2 



(7) 



which will be denoted by LJn. In particular, this poten- 
tial corresponds to a LJ(6,10) potential for n = 6. 

For HST potentials, it is not difficult to show that the 
K° parameter at zero energy is given by [0] 



Xhst(0,0 



^o(yo) cos(7U/ /2) - Y„ (y ) sin(7ri/ /2) 



Juo (Vo) sm(iris /2) + Y Vo (y ) cos(tti/ /2) ' 

(8) 

where v = (21 + l)/(n — 2), J and Y are the Bessel 
functions fTTJ, and y = [2/(n - 2)](/3„/n)) (n_2)/2 ■ The 
number of bound levels for angular momentum I is given 
by 



JVhst(0 = 



m ; 3vQ,m < UO < 3uo,m+l 
, VO < jv ,l 



(9) 



where j Vo ,m (jn > 1), is the m-th zero of the Bessel func- 
tion J Uo (x) @. 

For LJn potentials, the following results can be derived 
with the help of a local scaling transformation || . 

K£ Jn (0,0 = tam> /2)[l + hi(zo)][l - h^zo)]- 1 , (10) 



TABLE I: Selected data, all in atomic units, for effective po- 
tentials designed for the triplet state of a 23 Na dimer. All 
potentials have C 6 = 1556 a.u. Jl| and K c (0, 1 = 0) = 13.57. 
Ni = o is the number of s wave bound levels supported by the 
potential. D e is a derived parameter presented for discussion. 





LJ(6,10) 


HST 


A i=0 


Cio 


D e 


ro 


D e 


1 


1.65080e+9 


2.64709e-7 


3.22715e+l 


1.37752e-6 


2 


4.05415e+8 


2.17501e-6 


2.40025e+l 


8.13720e-6 


4 


1.00245e+8 


1.76896e-5 


1.74332e+l 


5.54309e-5 


16 


6.21105e+6 


1.14700e-3 


8.89948e+0 


3.13200e-3 



where z = (f3 n / f3 2n - 2 ) n - 2 /[2(n - 2)], 

= ^ o sin7r(z + 1/2 - ^o/2)r(z + 1/2 - i/ /2) 
l(Z0) Z ° sin7r(z + l/2 + ^o/2)r(zo + l/2 + j/ /2)' 

(11) 

and vo = (21 + l)/(n — 2). The number of bound levels 
for any I is given by 



N LJn (l) 



_ / [zo + \ - f ] , z > (yo + l)/2 
, z < (i/„ + l)/2 ' 



(12) 



where [x] means the greatest integer less than or equal to 
x. We note that scatte ring lengths, which are much more 
restrictive concepts [|[ |Il| , can be derived from K c (0, 1). 
For example, the s wave scattering length, which is de- 
fined only for n > 3, can be obtained from K c (0, 1 = 0): 



ai=o/p n 



.26 



r(i-6) 



r(i + 6) 



K c (0,0) + tan(7r6/2) 
K c (0,0) - tan(7ri>/2) 



(13) 



where b = 1/ (n — 2) . 

From these results, the HST and LJn types of po- 
tentials can be readily designed according to Eq. Q to 
have the desired K c (0,l), and the desired Ni. Table | 
gives a selected set of designs for the triplet state of a 
23 Na dimer. Here the potentials are designed to have 
K c (0, 1 = 0) = 13.57, which is found numerically us- 
ing the latest potential for a sodium dimer jl2). From 
Eq. ( p^| ) , this corresponds to an s wave scattering length 
of 64.57 a.u. The number of bound s wave levels sup- 
ported by this "real" potential is found numerically to 
be 16. 

Figure |l| shows the comparison of s and d wave par- 
tial cross sections of effective LJ(6,10) and effective HST 
potentials, both designed to support 16 s wave bound 
levels and have a K c (0, 1 = 0) = 13.57, with the Na-Na 
partial cross sections computed from the "real" potential 
p"2[ . The results are hardly distinguishable over a wide 
range of energies. In comparison, the hard-sphere poten- 
tial (HS) fails quickly away from the threshold for the s 
wave, and gives completely wrong results for the d wave. 

This result confirms the concept of effective potential 
based on AQDT. For our purposes here, however, what 
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- HST(16) 
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FIG. 1: Comparison of s and d wave partial cross sections 
of effective LJ(6,10) and HST potentials, both designed to 
support sixteen s wave bound levels and have K c (0,l = 0) = 
f3.57, with the "real" Na-Na partial cross sections E3] . The 
number in the parenthesis represents the number ors wave 
bound levels supported by a potential. Results for a hard- 
sphere (HS) potential are also shown for comparison. 
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FIG. 2: Comparison of s and d wave partial cross sections for 
an effective LJ(6,10) potential, designed to support a single s 
wave bound level, with the "real" Na-Na cross sections 



is more important is how shallow the effective potentials 
can be while still maintain a good description of low- 
energy characteristics of a real system. Figure || shows 
the comparison of s and d wave partial cross sections of 
an effective LJ(6,10) potential, in this case designed to 
support only a single s wave bound state, with the "real" 
Na-Na results j3. The agreements remain excellent. 

Figure |^ shows similar results for a HST potential. In 
this case, the HST potential that supports only a single 
s wave bound state does not do as well near the d wave 
shape resonance. But it quickly improves, monotonically, 
as the number of bound levels it is designed to support 
increases. By Ni = q = 4, a good agreement is achieved. 



FIG. 3: Comparison of s and d wave partial cross sections for 
effective HST potentials, designed to support one and four s 
wave bound levels, with the "real" Na-Na cross sections [|I4j. 
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FIG. 4: Zero-energy wave functions of effective potentials with 
four s wave bound levels compared with the "real" Na-Na 
wave functions flf4f. 



The robustness of these designs is not limited to the 
description of scattering properties, it also applies to the 
energies of bound states that are close to the dissociation 
threshold, and to the wave functions. For example, for 
effective potentials supporting four s wave bound levels, 
the HST gives a binding energy of 0.2027 GHz for the 
least-bound s state, the LJ(6,10) gives 0.2003 GHz. Both 
in good agreement with the result for the "real" potential, 
which gives 0.2044 GHz Q. Figure § shows that for 
effective potentials supporting four s wave bound levels, 
the wave functions are well represented down to r = 20 
a.u., covering basically all regions of space in which there 
is an appreciable amplitude. 

We stress that while only the results for sodium are 
presented here. They are used to illustrate much more 
general concepts. As the number of bound levels sup- 
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ported by an effective potential increases, all physical 
properties of states around the threshold converge to 
the same results (see Figs. [j]}4|). This is the shape- 
independence at length scale (3 n . The converged results, 
properly scaled, represent a set of universal properties 
shared by all quantum systems with the same type of 
long-range potential, and characterized by the same, 
/-independent constant K c = lim^/^^o K c (0, 1) 0, 
where (3 X represent the next shorter length scale present 
in the system. The examples presented here, Figs. ^ 
^, show how quickly this set of universal properties are 
approached as one increases the number of bound levels 
supported by an effective potential. This quick conver- 
gence is due to the fact that deviations from the uni- 
versal behavior depend on a high power of f3 x //3 n [||. 
Other quantum systems differ from Na primarily in K c 
JjJ , which does not effect this rapid convergence. 

Note that we did not make any distinction between 
two-atom and N-atom quantum systems in the state- 
ments above, because the same apply to a N-atom sys- 
tem. A short argument is simply that diffuse states, 
in which atoms are mostly at large distances relative to 
each other, only couple coherently to other diffuse states. 
A longer argument can proceed as follows. The correct 
K c (0,0), and therefore ai = o, ensures the correct results 
at the mean- field level ||. The correct two-atom wave 
function ensures the correct two-atom correlation. It 
also ensures the correct three-atom correlation, as fol- 
lows. Think of a three-atom as a two-atom perturbed by 
another. Frank-Condon considerations tell one that only 
two-atom states around the threshold are significantly 
coupled. This means by having the correct two-atom 
wave functions around the threshold, one has also the 
correct three-atom correlation around the threshold. . . . 

To illustrate the savings of computer resources as a re- 
sult of using an effective potential, consider the problem 
of N interacting atoms in a symmetric trap of frequency 
lu. If a real potential is used, the fact we need to rep- 
resent the length scale of (2fj,D e /h 2 )~ 1 ^ 2 means we need 
roughly D e /Tiijj number of harmonic oscillator states for 
each atom, for a total of (D e /h~u>) N number of states (ig- 
noring statistics). If an effective potential is used, the 
corresponding number is (D ete s/hui) N . Thus the sav- 
ing in the size of basis set is characterized by the factor 
(D e ,eff/ D e ) N . For the triplet state of Na, D ete g/D e is 
of the order of 10~ 3 if an effective potential with two s 
wave bound levels is used (see Table |) . This corresponds 
to a saving in the size of basis set of 10 9 fold just for a 
three-atom problem. Even greater savings are achieved 
for deeper potentials or for more atoms. From another 
angle, for effective potentials with A; =0 ~ 1, a ll length 
scales shorter than (3 n have effectively been eliminated. 
It is this elimination of short length scales that makes a 
complex problem more manageable. 

On the other hand, if a good description over an even 
wider range of energies around the threshold is desired, 



the same methodology can be carried to scales shorter 
than [3 n (e.g., /3$ for atoms with a —C^/r 6 long-range 
interaction and a — Cg/r 8 correction) B. However, be- 
cause the ratio Ps/Pe is different for different systems, the 
results become dependent upon one more system-specific 
parameter. At this stage, going to shorter length scale 
seems useful only in specialized two-atom applications 
We also point out that it is around the threshold 
that the quantum effects are most important |l5| . 

In conclusion, we have established the concept and de- 
sign of effective potentials describing atomic interaction 
at the length scale of j3 n . It is the scale that one has to 
deal with in studying quantum few-atom Q and quantum 
many-atom systems ||] at finite temperature, of high den- 
sity, or under strong confinement. We expect the method 
presented to play a role in our understanding of some 
of the more complex systems and processes at low tem- 
peratures, such as the three-body recombination process 
[I, excited clusters states ||, and quantum liquids 
[3|. In all cases, one can look for, and verify universal 
properties at the scale of /3„, by comparing results from 
different designs, such as HST and LJn, and by checking 
convergences as one relaxes an effective potential towards 
more bound levels. 
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